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The study of hot plasma expansion in a magnetic field is of interest for many laboratory and astrophysical 
applications. In this paper, an exactly solvable analytical model is proposed for an expanding resistive plasma 
in an external magnetic field in the regime in which the magnetic field does not perturb the plasma motion. 
The model is based on a class of exact solutions for the purely radial expansion of the plasma in the absence 
of a magnetic field. This approximation permits the reduction of the electromagnetic problem to consideration 
of a diffusion equation for the magnetic field. Explicit solutions are derived for a resistive cylindrical plasma 
expanding into a uniform ambient magnetic field. Some numerical examples related to the laser-produced 
plasma experiments are presented. 

PACS numbers: 52.30.-q, 03.50.De, 41.20.Gz, 52.65.Kj 



I. INTRODUCTION 

The problem of hot plasma expansion into a vacuum or into a background plasma in the presence of an external magnetic 
field has been discussed in the analysis of many astrophysical and laboratory applications (see, e.g., Refs. and references 
therein). In particular, such process is a topic of intense interest across a wide variety of disciplines, with applications to solar 
J6J and magnetospheric LI physics, astrophysics [9], and pellet injection for tokamak refueling ifioll . 

In this paper we consider analytically a resistive cylindrical plasma expanding into vacuum in the presence of an external 
magnetic field. Similar problem has been treated previously in Refs. lfTll[l2ll but for a spherical plasma expansion. Of the vast 
literature on the theory of plasma expanding into a vacuum or into a background plasma, Refs. 1 11-22] (see also references 
therein) illustrate various aspects and approaches. In Refs. ifTsUTsIl plasma has been considered as a highly conducting medium 
with zero magnetic field inside. From the point of view of electrodynamics, it is similar to the expansion of a superconductor in 
a magnetic field. An exact analytic solution for a unifornily expanding, highly conducting plasma sphere in an external uniform 
and constant magnetic field has been obtained in Ref. il3ll . The nonrelativistic limit of this theory has been used by Raizer 
lfl4ll to analyze the energy balance (energy emission and transformation) during plasma expansion. A similar problem has been 
considered in Ref. ifTsIl in a one-dimensional (ID) geometry for a plasma layer. In our recent papers, we obtained an exact 
analytic solution for the uniform relativistic expansion of a highly conducting plasma sphere fm Il7ll or cylinder [18^ in the 
presence of a dipole or homogeneous magnetic field, respectively. 

The mentioned treatments ifO - USll were obtained assuming a somewhat idealized situation: uniform expansion, infinite elec- 
trical conductivity of a plasma, etc. More realistic models for plasma expansion taking into account the deceleration (or accel- 
eration) of the plasma boundary have been developed, for instance, in Refs. 1 19-22] (see also references therein) for spherical 
ifigll . planar (ID) lf2(iil21i1 and cylindrical ll22ll expansions employing ideal magnetohydrodynamic (MHD) equations. However, 
it should be noted that the ideal MHD may not be justified in some experimental situations where the typical parameters are such 
that the plasma resistivity is not negligible 11,3, 11, 12]. In this case the coupling of the magnetic field with the plasma motion, 
determined by the magnetic Reynolds number, should result in a distortion and diffusion of the field across the expanding plasma 
lfTll [T2ll. We present here calculations of the electromagnetic field configuration in the stages preceding significant deceleration 
of the plasma and in the regime in which the magnetic field does not perturb the purely radial motion of the cylindrical plasma. 
The latter assumption is valid at large initial ratios of plasma energy density to magnetic field energy density. 



II. RESISTIVE MHD MODEL 



Usually, the motion of the expanding plasma boundary is approximated as the motion with constant velocity (uniform expan- 
sion). In the present study, a quantitative analysis of plasma dynamics is developed on the basis of a cylindrical model. Within 
the scope of this analysis, the nonuniform plasma expansion process is examined. We consider a resistive collision-dominated 
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magnetized plasma expanding into vacuum in the presence of a uniform and constant magnetic field. The relevant equations 
governing the expansion are those of resistive MHD [23J, assuming that the characteristic length scales for plasma flow are much 
larger than the Debye length and Larmor radius of the ions. Thus 



| + V.(pu) = 0, 

du . _^ 
p ^ + (u-V)u 



Vp = — [V X H] X H, 

47r 



— = V X [u X H] + DV^H, (1) 

with V • H = 0, where H is the magnetic field, D = c^/47rCT is the diffusion coefficient, p, u, -p and a are the mass density, 
the velocity, the pressure and the electrical conductivity of the plasma, respectively. In this paper we assume an isotropic 
and homogeneous electrical conductivity (and hence the diffusion coefficient D) of the plasma a. The equations above must be 
accompanied by the equation of state and the equation for entropy. Using the thermodynamic relation between entropy, pressure, 
and internal energy as well as Eq. ([T) the equation for pressure reads ll23ll 

^ + (u.V)p + 7P(V-u) = (7-l)^. (2) 
ot a 

Here 7 is the ratio of the specific heats^ is the current density in a plasma which after the elimination of some unimportant 
terms from the generalized Ohm's law ll23ll is reduced to the form 

j = a(^E+i[uxH]^ , (3) 

where E is the electric field. It is convenient to introduce the vector potential A. Within the scope of the present study the free 
charge density is absent and a suitable gauge V • A = allows the electric and magnetic fields to be determined from the vector 
potential A, 

E=— H = VxA. (4) 

c at 

Then from the last expression in Eq. ([T) one can derive a similar equation for the vector potential A 

(9A 

— = u X [V X A] + DV^A. (5) 

at 

In deriving Eq. (|5]l we have neglected the displacement current which is justified for the nonrelativistic expansion of the 
plasma. More specifically this approximation is valid at 47r(TTD ^ 1, where td = AnaR^ /(? is the characteristic diffusion 
time of the magnetic field and R is the characteristic size of the system, taken here as the radius of the cylindrical plasma. 
Alternatively, this inequality implies that cTt) is much larger than the plasma radius, ctd S> R, which is well justified for the 
nonrelativistic expansion velocities. 

The neglect of the Hall current in Eq. (|3) and the assumption of a scalar electrical conductivity for the magnetized plasma 
is justified when the characteristic time for the Coulomb collisions v^^ (where Vc is the collision frequency) is much less than 
the cyclotron period of the electron. Therefore, in this highly collisional regime the conductivity a is essentially a function of 
plasma temperature alone ll24ll . Hereafter it is assumed that the plasma temperature and the conductivity a{t) are uniform and 
are the functions only of time. 

The system of Eqs. ([rii-(|5]l can be used for determination of the dynamics of the expanding plasma as well as for the investi- 
gation of the evolution of the induced electromagnetic fields. For further simplification of this system we note that as long as the 
Lorentz force density ^ [j x H] and the Joule dissipation P /a terms are negligible compared with the hydrodynamic terms in the 
left hand side of Eqs. JTJ and dU, respectively, the plasma motion is a free radial expansion. For many realistic situations with 
laser-produced plasmas the free expansion is realized when the kinetic energy density of the expanding plasma is greater than 
the magnetic field energy density fr] (high-beta plasma), pv? /2 > Hq/Stt, where Hq is the strength of the initial unperturbed 
magnetic field. This is a necessary condition that the expansion remains radial and cylindrical. It has been shown [19| that the 
system of Eqs. ^ and (|2]i allows in this case the self-similar solutions for the quantities p, u, and p which are realized under 
specified initial conditions. These solutions are characterized by a radial velocity distribution linearly dependent on the radial 
coordinate r. At r ^ R{t), Ur{r, t) = r[R{t) / R{t)], where R{t) and R{t) are the radius and the velocity of the plasma boundary. 
In addition, the velocity Ur{r, t) vanishes at r > R{t), Ur{r, t) = 0. The self-similar solutions for the density p and the pressure 
p as well as the criterion of the violation of the free expansion solutions are discussed in Ref. [19]. However, we would like 
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to emphasize that since the hydrodynamic terms are reduced rapidly as the plasma expands and the Joule heating increases the 
plasma temperature and the electrical conductivity the electromagnetic terms in the right hand side of Eqs. ([1} and (|2) will no 
longer be negligible at the final stage of the plasma expansion when plasma may fully be stopped and deformed by the magnetic 
field pressure. As mentioned above, the average plasma pressure p is strongly reduced compared to the magnetic pressure and 
the model of the purely radial expansion clearly becomes invalid in this case. Nevertheless, if the critical time interval At, where 
p < H'^ /Stt, is much smaller than the typical time scale of the plasma flow (up to the full stop), the contribution of this interval 
to the overall plasma dynamics is negligible and use of the radial expansion model is justified. 

In the next sections we will use the profile of the plasma radial flow velocity Ur{r, t) — r[R{t) / R{t)] together with Eq. (|5]l to 
investigate the electromagnetic field configuration generated by the expanding cylindrical plasma. 



III. SOLUTION OF THE MOVING BOUNDARY AND INITIAL VALUE PROBLEM 



In this section we consider the moving boundary problem of the plasma cylinder expansion in the vacuum in the presence of 
the constant and homogeneous magnetic field Hp. Consider a cylindrical region of space with radius r = R{t) at the time t 
containing a neutral plasma which has expanded at t = (with R{0) = Rq) to its present state from a cylindrical source with 
radius Rq located around r = 0. We assume that at any time t the plasma cylinder is unbounded in z direction (i.e. the plasma 
cylinder is located at — oo < z < oo). To solve the boundary problem we introduce the cylindrical coordinate system (r, if, 
z) with the z-axis along the plasma cylinder symmetry axis and the azimuthal angle ip is counted from the plane (a;z-plane) 
containing the vector of the unperturbed magnetic field Hq. The angle 6 between the vector Hq and the z-axis is arbitrary. 

As the cylindrical plasma expands it both perturbs the external magnetic field and generates an electric field. We shall obtain 
an analytic solution of the electromagnetic field configuration. We consider the case of the purely radial expansion of the 
plasma cylinder with an arbitrary (but nonrelativistic) expansion velocity R{t). Having in mind the symmetry of the unperturbed 
magnetic field and the fact that the electromagnetic fields do not depend on the coordinate z it is sufficient to choose the vector 
potential in the form Ar — 0, 

Ho\\W{r,t) = ^M^), A, = i/o±*(r,Osin^, (6) 

where W{r, t) and ^(r, t) are some unknown functions. From symmetry considerations the functions W{r, t) and ^'(r, t) are 
independent on the cylindrical coordinate tp. Here i7o± and iJo|| are the components of the unperturbed magnetic field Ho 
transverse and parallel to the z-axis, respectively. The components of the electromagnetic field are expressed by these functions 

as if, = Foil 



— cos ^p, ^ -Ho±_ — 

r or 



Hr = Ho± — COS (fi, = -Ho± — sin ip, (7) 



1 dA^ 1 9* 

E^^ ^, E, = — i7o^_ sinvj, (8) 

c ot c ot 

and Er = 0. The equation for the vector potential A(r, t) inside the plasma cylinder (r ^ R{t)) is obtained from the MHD 
equation for the magnetic field diffusion, Eq. Q, which for the unknown functions 'I'(r, and W{r,t) yields a system of 
equations 

dt ^ dr \ dr'^ r dr / 

^j^^i^uW)- d(^^--— + —\ (10) 
dt dr ^ \ dr'^ r dr r^ J 

Here Ur{r, t) = r[R{t) / R{t)]. In the vacuum surrounding the plasma cylinder (r ^ R{t)), the magnetic field is determined by 
the Maxwell equation V x H = with j = which for the functions '^{r, t) and W{r, t) become 

15* * 0, (II) 



Qj,2 J. Qj, y,2 

d^W 1 dW W 



^=0. (12) 



Qj,2 J, Qj. J, 

In the plasma the magnetic field is the solution of the diffusion equation with a diffusion coefficient D{t) = c^/ATTa{t). The 
plasma conductivity a{t) is essentially a function of plasma temperature alone. If it is assumed that the plasma temperature is 
uniform, this coefficient is a function only of time and Eqs. (l9]l-(fT2]i may be solved by the method of separation of variables. 
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The system of equations (l9]l-(fT2]i is to be solved in the internal (r ^ R{t)) and external (r ^ R{t)) regions subject to the 
boundary and initial conditions. Since the plasma under consideration has finite electrical conductivity, there are no surface 
currents at the plasma-vacuum boundary and H must be continuous at the moving boundary. The continuity of the magnetic 
field H requires that 'i/{r,t),W{r,t) and the radial derivative d'^/drhe continuous at the expanding plasma surface. In addition 
we require that as r — >^ oo, the magnetic field H(r, (p, t) be time-independent and asymptotically approach the uniform magnetic 
field Ho. This is equivalent to the boundary conditions ^{r,t) = W{r,t) = r at r — > oo. 

The initial conditions are at < = 0. In principle two distinct sets of initial conditions could be considered ifTTl] . 

(i) In the case of the poorly conducting plasma the initial conditions are imposed for arbitrary r: 

^t{r,0) = Wir,0)=r. (13) 

(ii) In the case of the perfectly conducting plasma the initial conditions are imposed separately for the domains inside (r ^ Rq 
with Rq — R{Q)) and outside (r ^ Rq) the plasma cylinder [18]: 

^ {r,0) = W {r,0) = 0, r sC i?o, 

^(r,0) = r ^, W(r,0)^r, r ^ Rq. (14) 

r 

In the first case, Eq. (flJl i. the plasma is initially called and poorly conducting so that the external magnetic field is completely 
penetrated inside the plasma. At the other extreme case (ii), the initial conditions (fT4l) imply that initially the plasma is highly 
conducting so that the magnetic field is completely excluded from the initial plasma volume. We consider below the initial value 
problem (ii). The extension of the obtained solution to the case (i) is straightforward. 

At r ^ R{t) we look for the solutions of Eqs. (fTTT i and (fT2] i for the functions ^'(r, t) and W{r, t) in the form ^ r", where 
a is some numerical constant. Therefore, taking into account the boundary condition at r — ^ oo the full solution in the domain 
outside the plasma cylinder is given by 

^>{r,t)^r~C{t)^, (15) 



r 



W{r,t) = r 



1 + Ci(t) In 



r 



R{t) 



(16) 



where C{t) and Ci{t) are the arbitrary functions of time with the initial conditions C(0) = 1 and Ci(0) — 0. However, since 
the magnetic field should be finite at r — > oo we set Ci{t) = 0. At r ^ R{t) this gives the final solution W{r, t) = r for the 
function W{r, t). 

For the class of separable solutions the motion is such that, for a given element of plasma, the quantity ^ = r/R{t) is a 
constant of the motion. The solution of Eqs. Q and ( fTol ) inside the plasma cylinder (i.e. at r ^ R{t)) is then facilitated by the 
representation of the functions 'I'(r, t) and W{r^ t) in the form 

oo 

^'(r,i) = r + ^a„r„(i)$„(C), (17) 

^ oo 

W{r,t) = r + —-Y,bnUn{t)Qn{0. (18) 

^ ^ n—1 

where T„(i), Un{t), and 8„(^) are some unknown functions, and a„ and 6„ are the unknown expansion coefficients. 

Next inserting these expansions into Eqs. (|9]l and ( fTOl i for "I>n(C) ™d one arrives at the ordinary differential equations for 

the cylindrical functions ll25ll . We choose only the regular solutions of the obtained equations finite at the origin r = (or at 
^ = 0). Thus 

Here An and Bn are the integration constants, A„ and k„ are some arbitrary parameters (depending only on n), arising due to 
the separation of the variables, and Jq and Ji are the Bessel functions of the first kind. 

In the same way for the time-dependent functions Tn{t) and Un{t) we obtain the following set of the ordinary differential 
equations: 

tn{t) + XlmTn{t)^R{t), (20) 
Unit) + nlmUnit) = 2R{t)R{t), (21) 
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where 



D{r) 



dr. 



The general solutions of the first order differential equations (|20] | and (ISTT l can be represented in the form 



(22) 



Tn{t) 



Unit) 



to + I e 





R{T)dT 



uq + 2 I e 





R{T)R{T)dT 



(23) 
(24) 



Here to — Tn (0) and uq — [/„ (0) are the initial values of T„ (t) and J7„ (t), respectively, to be determined by imposing the initial 
conditions. In addition, it should be emphasized that the expansions given by Eqs. ([TtT i and ( fTsT i are the solutions of the system 
of Eqs. (|9]l and (fTOl i only if the functions and QniO satisfy at arbitrary ^ ^ ^ 1 the equations 



n—1 n—1 



(0 - 



(25) 



These relations impose some additional constrains on the expansion coefficients a„ and 6„. 

Now let us consider the boundary condition at r = R(t) for the functions W{r, t), '^(r, t), and ^^'(r, t). For the function 
W{r, t) this boundary condition yields Jo{Kn) — 0, i.e. the quantities k„ (with n = 1,2, ...) must be the positiv zeros of the 
Bessel function Jo{z). Later on we will assume that the zeros k„ are arranged in ascending order of magnitude. The same 
boundary condition for the quantity 5'(r, t) determines the unknown function C{t) in Eq. (fTsl l. 



Cit) 



^ n—1 



a„A„r„(t) Ji(A„). 



(26) 



Finally, the boundary condition at r = R{t) for -g:p'^{r, t) yields another relation for the function C{t) which should be consistent 
with Eq. (l26T l. This is only possible if Jo(A„) — 0, i.e. A„ = k„ are also the zeros of the Bessel function. 

Let us now turn to the determination of the expansion coefficients a„ and 6„ using the constrains in Eq. ( l25l l. Inserting Eq. (fT9l l 
into Eq. (l25l l and using the summation formulas (IA5b and (IA2b one can easily prove that the equations in Eq. (l25l l are satisfied 
if a„A„ = — 4/[A^j Ji(A„)] and 6„i?„ = — 2/[A„ Ji(A„)]. Then taking into account the summation formulas ( IA6b and ( IA2b 
the initial conditions C(0) = 1 and W{r, 0) = (at r ^ i?o) for the functions C{t) and W^(r, i) imply that T„(0) = i?o and 
Un{^) = Rq^ respectively. Moreover, having in mind the relations (|A5) and (|A6b the initial conditions for the quantity 5'(r, 0) 
inside (r ^ Rq) and outside (r ^ Rq) the plasma cylinder (see Eq. (fl4] i) are then satisfied automatically. 

Therefore, the complete solution of Eqs. (l9ll-(fT2t subject to the initial and boundary conditions inside (r R{t)) and outside 
(r ^ i?(i)) the plasma cylinder, respectively, may be represented as 



n=l 



Ji(A„0 
A2Ji(A„)' 



W{r,t) = r 



1 —\^U fi)jMML 



and W{r, t) = r. 



(27) 
(28) 

(29) 



with 



oo ^ oo 



A2 



(30) 



Note that \l'(r, i) — W{r, t) = r at r — > oo as expected by the boundary conditions at the infinity. 
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The (^-component of the vector potential is determined by the first relation in Eq. (|6]l. The straightforward integrations in 
Eq. ( |28] ) with respect to the radial coordinate r result in ||2^ 



G{t) 



(31) 
(32) 



inside (r ^ R{t)) and outside (r ^ R{t)) the plasma cylinder, respectively. Here G{t) is an arbitrary function of time to 
be determined by the boundary condition at r = R{t). From this condition one obtains G{t) — 4U{t). Equations (|27|i-(|32]i 
represent the complete solution of the problem and determine the structure of the electromagnetic fields both inside and outside 
the expanding plasma cylinder. Expressions for the components of the electromagnetic fields E and H may now be obtained by 
use of Eqs. (|7]l and (|8]l. These components are, for r ^ R{t), 



Hr = Hq± cos (p 

= -Hoi_ sin (fi 
Hz = Hq\\ 1 
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A2Ji(A„)_ 



mt{ A„ji(A„) 



(33) 
(34) 
(35) 



and, for r ^ R{t), = Hqu, 



2H. 



oil 



CO ^ 



Un{t)MKO 



m 
m 



H{t) 



Hr — Hq± cos if 



1 - «r(t) 



J,. 

AR{t) 



m 



(36) 



(37) 



(38) 
(39) 



o|| 
cr 



R{t)R{t) 



2D{t) 



sinip{R{t)[R{t)+A:T{t)] 



Am 

R{t) ^ 

^ ^ n—1 



(40) 
(41) 



In the latter expressions for the components of the electric field the time-derivatives Un{t) and Tn{t) have been excluded by 
means of Eqs. (|20] | and (ISTT i. Also in Eqs. (l33Tl-(l4ni the prime indicates the derivative of the Bessel function with respect to the 
argument. 

The induced current density j is defined only inside the plasma cylinder (i.e., for r ^ R{t)) by the relation (|3]l (or using the 
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Maxwell equation j 



4lT 



^ V X H) and has the following components: jV = 0, and 



oo 



Jl(AnO 
Jl(A„) ' 



27ri?3(t) 



(42) 



n=l 



7ri?2(t) 



oo 



Ji(A„) ■ 



(43) 



n=l 



Until now we have considered the initial value problem (ii) assuming that initially the plasma conductivity is so high that the 
magnetic field is completely excluded from the initial volume of a plasma. The extension of the obtained solution to the case of 
the initial value problem (i) (with highly resistive plasma at t = 0) is straightforward. From the consideration above it follows 
that the solution of the boundary and initial value problem (i) is again determined by Eqs. (|27]|-(|43T|. where, however, the initial 
conditions T„(0) = Rq and Un{0) = Rq for the functions Tn{t) and Un{t) (see Eqs. (l23T l and ( l24b ) should be replaced by the 
zero initial conditions, T„(0) = Un{0) = 0. Comparison of the complete solutions obtained by the initial value problems (i) and 
(ii) shows that the electromagnetic fields and the induced current for the two distinct cases differ only in the terms containing 
r„(0) and Un{0) in Eqs. ( |23] | and ( l24l i. These terms force the matching of the solution to the initial condition of the magnetic 
field completely excluded initially from the plasma volume r ^ Rq. Because of their exponential dependence on d{t) (and 
hence on the time t) these terms become negligible compared to the other terms in the electromagnetic fields in the time required 
for 'd{t) to become of the order ■d{t) ^ 1. This time interval is the characteristic time for diffusion of the magnetic field into a 
stationary plasma of radius R and conductivity cr, i.e., td — R^ /D = AiTaR^/c^. Thus the initial conditions for the initially 
perfectly conducting plasma are forgotten by the plasma at < > tq. 

In the context of the two distinct initial value problems (i) and (ii) it should be also emphasized that the initial value of the 
plasma conductivity cr(0) should be consistent with the chosen physical model. Indeed, the cases (i) and (ii) imply vanishing 
((t(0) — > 0) and very large ((7(0) — ?> oo) initial conductivities of the expanding plasma, respectively. As a demonstration of the 
importance of the initial value a{0) consider, for instance, Eqs. (|36] |. (l37t and (l40l i. (HTt for the generated electric field. Using 
Eqs. ( I20I 1 and ( 1211 1 as well as the summation formulas of AppendixlAlit is straightforward to show that at t = the electric field 
inside (Eqs. ( |36] | and ( |37] |) the plasma cylinder is given by 



where Rq ~ R{0). Thus, the initial electric field in the plasma volume vanishes or is finite in the cases (ii) (with uq = Rq, 
to = Rq) and (i) (with uq ~ to ~ 0), respectively. The initial electric field in a vacuum is determined by Eqs. ( |40l i and ( 1411 1 at 
t = 0. It is seen that at t ^ the last terms in these expressions proportional to the diffusion coefficient D{t) vanish for the 
initial condition (i) while diverging as ^ in the case (ii). In the latter case assuming, however, a perfectly conducting 

(in fact infinitely conducting) plasma at t = one must consider the limit D{0) — > in the last terms of Eqs. ( l40l i and (I4l1 i 
which vanish eventually at t — > 0. Finally, we note that at f — > the nonvanishing terms in Eqs. ( |40l i. (l4T1 i and (l44l l. ( |45l l 
are proportional to the initial expansion velocity Rq of the plasma. Therefore, it is not surprising that at Rq ^ the plasma 
expansion builds up instantly an initial electric field although the induced magnetic field is zero. 

At the end of this section consider briefly the nondiffusive limit of the obtained solutions, Eqs. (I33ll-(l43ll. when the diffusion 
coefficient vanishes, D ^ Q. This limit can be obtained using at 13 ^ the expressions (I22li-(l24li which yield r„(<) — R{t) 
and Unit) = R^{t). Therefore, having in mind the summation formula (|A6t from Eq. ( [30l l one finds T{t) = R{t)/A and 
U{t) — R^{t)/A. Using these results and the summation formulas derived in Appendix lAl it is straightforward to show that 
jip — iz = ^ and the electromagnetic fields are, for r ^ R{t), H(r,t) = E(r,t) — 0, and, for r ^ R{t), = i?o||' 
Hr = i/o_L cos (^(1 -i?Vr2), = -Ho±sh-np{l + R"^ /r^), ^ l3Ho\\{R/r), = 2pHo±{R/r) simp, where /? = R/c. 
These expressions have been derived previously in Ref. fil8l1 for the cylindrical plasma expansion neglecting the diffusion of the 
magnetic field. 




(44) 




(45) 



IV. ENERGY BALANCE 



Previously significant attention has been paid 10, IH [Til - fisll to the question of what fraction of energy is emitted and lost in 
the form of electromagnetic pulse propagating outward of the expanding plasma. In this section we consider the energy balance 
during the plasma cylinder expansion in the presence of the homogeneous magnetic field. 
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(46) 



Our starting point is the energy balance equation 

where S = ^ [E x H] is the Poynting vector and j is the induced current. Note that the density of the electric field energy has 

been neglected in Eq ( l46l l since R c and E <^ H. The energy emitted to infinity is measured as a Poynting vector integrated 
over time and over the lateral surface Sc of the cylinder with radius Vc, length Ic and the volume fic (control cylinder) enclosing 
the plasma cylinder (rc > R{t)). Integrating over time and over the volume ilc Eq. ( |46] | can be represented as 

Ws{t)^Wj{t)+AWM{t), (47) 



where 



Ws{t)=r, / dt' / Srd^, 
Jo Jo 

W]{t)^-^ f dt' f j-Edr. 



(48) 



Here Sr is the radial component of the Poynting vector Note that the total flux of the energy over the bases of the control cylinder 
determined by Sz vanishes due to the symmetry reason. Wm {t) and AT4^m (t) — Wm (0) — Wm {t) are the total magnetic energy 
and its change (with minus sign) per unit length in a volume ilc, respectively. Wj{t) is the energy (per unit length) transferred 
from plasma cylinder to the magnetic field and is the mechanical work with minus sign performed by the plasma on the external 
magnetic pressure. At i = the magnetic fields are given by H = Ho in the model (i) and by Eq. ( fT4l ) in the model (ii). Hence 



in (i) Wm (0) is the total magnetic energy per unit length in a volume flc and is given by Wm (0) = Qc 



r^(H^/87r) while in 



(ii) Wm (0) is the total magnetic energy in a volume Q'^ of the control cylinder excluding the volume of the plasma cylinder (we 
take into account that H = 0at<: = 0ina plasma cylinder in the model (ii)) and Wm(0) — Qd^ — cos^ 9 — k^** sin^ 9). 
Here k — rd Ro and is the angle between Hq and the symmetry axis of the plasma cylinder (z-axis). Then the change of the 
magnetic energy AWM(i) in a volume Q,c can be evaluated as 



AWmW 



(49) 



Hence the total energy flux Ws(i) given by Eq. ( 1481 1 is calculated as a sum of the energy loss by plasma due to the external 
magnetic pressure and the decrease of the magnetic energy in a control volume ilc- For expansion of a one-dimensional plasma 
slab and for uniform external magnetic field Ws — 2 Wj ~ 2AWm, i.e., approximately the half of the outgoing energy is gained 
from the plasma, while the other half is gained from the magnetic energy ifTsll . In the case of expansion of highly conducting 
spherical plasma with radius R in the uniform magnetic field Ho the outgoing energy Ws is distributed between Wj and AWm 
according to Wj — l.SQo and AWm = 0.5Qo with Ws = ^Qq, where Qo — H^R^ /Q is the magnetic energy escaped from 
the spherical plasma volume iflill . Therefore, in this case the released magnetic energy is mainly gained from the plasma. 
Consider now each energy component Ws(t), Wj(t) and AWM(i) separately. Ws(t) is calculated from Eq. (|48l l using the 



expressions for the electromagnetic fields generated outside the plasma (r ^ R{t)). RecaUing that Hz 
Eqs. (|39]l-(|4T]i and (|48]l we obtain 



Hqu at r ^ R{t) from 



Wsit) - i?o'|| 



2R^t) 



to 



o± 



(50) 



Here T{t) andU{t) are given by Eq. dSOl l. 

Next, we evaluate the energy loss Wj (t) by the plasma which is determined by the induced current density, j. This current has 
two azimuthal and axial components and has been determined in Sec. |III1 see Eqs. (|42] | and ( |43] l. Since the current is localized 
only within a plasma volume flu, in Eq. (l48T l the volume ilc can be replaced by the plasma volume JIr. The total energy loss by 
the plasma cylinder is calculated as 



Wj{t) 



m 

R^t) 



ARl 



+21 

''0 



(51) 



where 



lit) 



n=l " 



T^{t), m 



n=l " 



urn- 



(52) 
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FIG. 1: (Color online) The ratios Ts{t) (solid lines) and rj(t) (dashed lines) for k — 10 and 5 = 0.1 as a function of t (in units of 
tr — Ro/v) calculated from expressions ( 150) and ( 151) with ^ = (left panel), 6 — n/i (middle panel), and 6 — 71/2 (right panel). The 
dotted lines represent the quantity FiviCi) = rs(t) — rj(t). 



The change of the magnetic energy in a control cylinder is calculated from Eq. ( |49] l. For evaluation of the magnetic energy 
inside and outside the plasma volume we use Eqs. (|33]|-(|35]| and (l38l l, (|39] | together with Hz = Hq^\, respectively. Thus the 
change of the total magnetic energy in the control cylinder is represented as 



jj2 



2U{t) - 



m 



Up ( ^ Up 

4 Rl 



+Hl^lR{t)T{t)-%{t)+2T\t) 



R^jt) _ t_ 
R^K-^ i 



+ 2 {Rq — t(j 



(53) 



Comparing now Eqs. (l50t , (ISTT i and ( l53T l we conclude that AVFM(i) + W,]{t) = W^{t) as predicted by the energy balance 
equation ( |47] i. Let us recall that Eqs. dSOl l, ( ISTT i and ( |53T l are valid for both initial conditions (i) and (ii) choosing appropriate 
values for the quantities and uq (i.e., = uq ^ i) in (i) and = Rpi = (ii))- 

The energy components W^s it) and A Wm {t) depend naturally on the radius Vc of the control cylinder while Wj [t) is deter- 
mined only by the volume of the plasma cylinder. At rc —> oo Eqs. dSOl l and ( |53] | are finite and in this case Ws{t) represents 
the electromagnetic energy emitted to infinity. It should be noted that at r cx) the induced magnetic and electric fields behave 
as H ^ r^^ and E ^ r^^, respectively, see Eqs. dSSll- dTTl l. Consequently, at rc —> oo the induced magnetic field does not 
contribute to the energy flux VFs [t) and the latter is determined by the electric field and the unperturbed magnetic field Hq. 



V. NUMERICAL EXAMPLES 



In this section using theoretical findings of the preceding sections, we present the results of our model calculations for the 
electromagnetic fields generated due to the radial expansion of the magnetized cylindrical plasma into a vacuum. As an example 
we consider the radial expansion with R^ (t) = R^+v^t^, where v is the expansion velocity ait ^ co. Note that initially Rq = Q 
in this case. The electrical conductivity is modeled by a{t) = ap{T/t), where t is some characteristic decay time of (j{t) and 
(To is a constant. Accordingly the diffusion coefficient is given by D{t) = DqII/t) with Dq = c^/Anap. Thus, initially the 
conductivity of the plasma is very high which corresponds to the initial condition (ii). However, as plasma expands the plasma 
will be eventually cooled off and the conductivity decreases with time. For the chosen model 'd{t) — 5 \n['ri{t)] as it follows from 
Eq. ( |22] |. Here 6 — Dq/v^t is the dimensionless diffusion coefficient while ri{t) — R{t)/Ro is the dimensionless plasma radius. 
In order to obtain a physical idea of the length and time scales involved, let us consider briefly a numerical example. Taking, for 
instance, Rq = 1 mm and v = 10^ cm/s, one obtains tr = Rp/v = 0.1 ^s. For a laser generated plasma with a temperature 
T ~ 10^ K, one can take ctq — 5 x 10^** s^^ and Dq ~ 2 x 10^ cm^/s. Using these parameters the dimensionless diffusion 
coefficient S ^ 1 implies the characteristic decay time r ~ 0.18 /is for the conductivity. 

Equations (l23T l and (l24l i can be evaluated in a closed form 

"^^"R^ ° — KsTl — ' ^^^^ 

"^^^ra^ ° ■ ^^^^ 
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FIG. 3: (Color online) The strength of the magnetic field (in units of Ho) at r = as a function of time t (in units of tr) for 5 = 0.1 (left 
panel), 5 — 1 (middle panel), and 5 = 5 (right panel). The solid, dashed and dotted lines correspond to 6 = 0, 8 = n / 4 and 8 = n/2, 
respectively. 



It is seen that the first terms of Eqs. (l54t and (ISST i containing the initial values to and uq vanish exponentially at t ^ oo. The 
electromagnetic fields, the induced currents and the components of the energy are determined by inserting Eqs. (l54l and (l55l l 
into Eqs. (l33Tl-(l43]l and Eqs. (I50ll-(l53]l. respectively. 

First, we consider the energy balance during the plasma expansion. As an example in Figs.[T]and|2]we show the results of our 
model calculations for the ratios Tsit) = Wsit)/Q{t), T]{t) = Wj{t)/Q{t) and Tuit) = AWM{t)/Q{t) = Tsit) - rj(t) 
as a function of time (scaled by tr = Ra/v) for k = rc/Ra = 10, and for various orientations of the magnetic field Hq, 
9 = 0,9 = 7r/4, and 6 = tt/2. Let us recall that Eqs. ^ and ^ are vaHd at ^ t/rR ^ (k^ - l)i/2 (or r]{t) ^ k). Here 
Q{t) ~ {ttR^ (t) ) {Hq /Stt) is the total energy (per unit length) of the unperturbed magnetic field at time t in the volume of the 
plasma cylinder For the dimensionless diffusion coefficient 6 we have chosen two distinct values, (5 = 0.1 (weak diffusion) 
and 6 = 1 (strong diffusion). Unlike the cases with uniform magnetic field and highly conducting spherical or one-dimensional 
plasmas mentioned in Sec.|IV](see also Refs. ifTilfTsIl ) there are no, in general, simple relations between the energy components 
Ws{t), Wj(t) and Q{t). Moreover, in the present case with a finite conductivity of the plasma the ratio rj(t) is not a constant 
anymore (see Figs. [T]and|2]i as predicted in Ref. |18] for highly conducting cylindrical plasma. It has been shown in Ref. |18] that 
at cr oo the energy flux Ws (<) is systematically positive and F j (<) ^ 1, i.e. the energy loss by a plasma is approximately equal 
to the magnetic energy escaped from the highly conducting plasma volume (formation of a diamagnetic cavity). As demonstrated 
in Figs. [T] and |2]both energies Wj and Ws are negative here at the initial stage (at t < tr) of the plasma expansion. That is, 
the flux of the electromagnetic energy over the surface of the control cylinder is negative and the plasma gains energy from the 
generated electromagnetic fields. It is notewoithy that the energy gained by the plasma is maximal for transverse orientation of 
Hq. The physical origin of these features is the diffusion of the magnetic field in an expanding plasma. Alt = the magnetic 
field is completely excluded from the plasma volume. However, as soon as the plasma expands and the electrical conductivity 
is reduced the magnetic field is effectively diffused into the volume of the expanding plasma giving rise to the negative energy 
flux over the control surface. Clearly this process is intensified by increasing the diffusion coefficient as shown in Fig.|2] At 
the later time, t > tr, with essentially reduced a{t) the magnetic field almost freely fills the whole volume of the plasma and 
both energies Wj and Ws become positive in this case. That is, the expanding plasma loses its energy due to the transfer of the 
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FIG. 4: (Color online) The distribution of the magnetic field strength (in units of Ho) in the xy plane att = O.Itr for 5 = 1 and for S = 
(left panel), 9 = n/4 (middle panel), and 9 = n/2 (right panel). The coordinates are scaled by _Ro. 




x/R^ x/R„ x/R, 

FIG. 5: (Color online) Same as in Fig.|4]but att — tr. 



energy to the electromagnetic fields. At i > tr the quantities Tj{t) and Ts{t) are determined explicitly by 

r J (t) ~ Sa {6/2) cos2 e + 2S2 (5) sin^ 6*, (56) 
Ts (t) ~ 2Si (,5/2) cos^ 61 + S3 (S) sin^ 0, (57) 

where S3(^) = 25i(5) + (?7^/k^)S^((5) and the functions and 22(2) are given by Eqs. jAlOl l and (lAl 11 1, respectively. 

Note that the energies Wj{t) and Ws{t) are reduced by the magnetic field diffusion which at 5 > 1 results in Wj{t) ^ Ws{t). 

Next, we consider the temporal and spatial distributions of the magnetic field for various diffusion coefficients S and for 
orientation angles 9. Figure |3] demonstrates the absolute values of the magnetic field H{r, t) (in units of Hq) on the symmetry 
axis of the plasma cylinder, r = 0, as a function of time. As it follows from Eqs. ([33]l-(l37]i the electric field vanishes on this 
axis, E = 0, while H,. = Ho±Ti{t) cos (p, = —HQ±Ti{t) sin (p, = i?o||-^2(i), where 

Thus, as expected, the magnetic field H{r, t) is independent on the angle ip. From Fig.|3]it is seen that simultaneously with the 
plasma expansion the magnetic field inside grows monotonically from zero value and saturates typically at tsat — Rq/ Dq ^ S^^. 
Clearly the saturation time of the magnetic field decreases with S. At i > tgat the magnetic field is constant because it is 
completely diffused into the expanding plasma remaining, however, smaller than Hq. In this asymptotic regime from Eqs. (|54] i 
and (|55]) one obtains T„(t) ~ R{t)/{X'^S + 1) and Unit) ~ R^{t)/{\\5/2 + 1). Insertion of these expressions into Eqs. (|58] l 
and ( |59] l determines the magnetic field at t > ts^t ■ It is also noteworthy the influence of the orientation of Hq on the magnetic 
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FIG. 6: (Color online) Same as in Fig.|4]but at t = 2rR. 



field inside the plasma cylinder at r = 0. At weak diffusion, (5^1, the quantities JFi{t) and T2{t) are evaluated using Eq. ( |A7| i 
with a; = and z = S-^/'\ In this case we obtain J^i (t) ~ {Att^ /Sy/'^ exp{-l/ VS) and J"2(0 =i {Sn"^ /S)^^^ exp{- It 
is seen that J^i(t) ^ -7^2(0 and the magnetic field is essentially larger for the transverse orientation of Hq (Fig. [3] left panel). 
However, in the strongly diffusive regime, S ^ 1, from Eqs. ( |58] l and ( |59] l one derives J^i (i) ~ 1 — 1 /45 and T2{t) ~ 1 — 1/ 2(5, 
and the magnetic field is only weakly sensitive to 9 approaching the value Ho of the unperturbed magnetic field (Fig. |3] right 
panel). 

More complete spatio-temporal evolution of the magnetic field is shown in Figs. |4}{6] Shown are the distributions of the 
magnetic field strength in the a;y-plane transverse to the symmetry axis of the plasma cylinder (z-axis) at t — O.Itr (Fig.|4|i, 
t = tr (Fig. |5]i, and t ~ 2rR (Fig. |6]l and for 6 = 1, 6 = 0, 9 — n/A, and 6 = tt/2. Note that the plasma radii are 
R{t) ~ 1.005i?o, R{t) ^ lARo, and R{t) ~ 2.2i?o in Figs. |4] |5] and|6] respectively. It is seen that the magnetic field 
distribution essentially depends on the orientation angle d of the unperturbed magnetic field Hq. At the initial stage of the 
plasma expansion the magnetic field inside is excluded for all 6 except the domains in the close vicinity of the plasma boundary 
(r ~ R{t)) where the magnetic field diffuses into the plasma due to the finite conductivity. Moreover, outside the expanding 
plasma and at parallel orientation with 9 = the magnetic field is not perturbed (H = Hq) while at the oblique (6 = 7r/4) and 
transverse (6 = tt/2) orientations the magnetic field distribution exhibits much more complex behavior It is clearly seen that 
the expanding plasma pushes the magnetic field out in the direction of the transverse component Hq^ (i.e. at 1^9 = or = tt) 
where the strength of the magnetic field exceeds essentially the unperturbed value Hq. This effect is maximal at9 = tt/2 when 
H ~ 1.9Hq (see Fig.|4] right panel) in the vicinity of the plasma boundary at = and ip = tt. In the plane transverse to Ho± 
and containing the z-axis (i.e. at ip = 7r/2 and ip = 3tt/2) the expanding plasma forms a magnetic cavity outside its volume 
where the strength of the magnetic field is reduced and H ~ O.SHq (see Fig. H] middle and right panels). Thus, at this early 
stage of the expansion of the highly conducting plasma the role of the magnetic field diffusion is small. 

Next, at the later time, t > tr, the electrical conductivity of the plasma is reduced and the magnetic field is essentially 
diffused into the plasma which is more effective at larger 6 (Figs. |5] and |6]l. Furthermore, in the domain outside the expanding 
plasma volume the diffusion smooths out the magnetic field distribution and it is only weakly perturbed for parallel and oblique 
orientations of Hq. At transverse orientation, however, the magnetic field deviates considerably from Hq. In particular, it is 
noteworthy that in the latter case the magnetic cavity is compressed along the direction of Ho± (y-axis). Finally, att> 2tr 
with a further reduction of the plasma conductivity the magnetic field diffuses considerably into the plasma (Fig.|6]l eventually 
recovering the uniform distribution in the xy-plane. However, the strength of the magnetic field inside the plasma volume is 
systematically smaller than Hq as shown in Fig. [3] 



VI. CONCLUSION 



In this paper, an exact solution of the purely radial expansion of a neutral, resistive plasma cylinder in the presence of a 
magnetic field has been obtained. The electromagnetic fields are derived by using the appropriate initial and boundary conditions. 
Two kind of the initial conditions have been considered assuming poorly or perfectly conducting plasma at the initial state. In 
the first case the external magnetic field is completely penetrated inside the plasma while in the other extreme case the magnetic 
field is completely excluded from the initial plasma volume. However, as expected both solutions "forget" the imposed initial 
conditions at t > td, where td is the characteristic diffusion time. Using a simple model for the electrical conductivity we have 
also studied the energy balance during the plasma expansion as well as the spatio-temporal distribution of the magnetic field 
with an arbitrary orientation of the initial field. In contrast to the previous treatments with highly conducting plasmas ifTsUlsll . 
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our model calculations demonstrate some new features arising due to the finite conductivity of the expanding plasma. 

Going beyond the present model, which is based on several approximations, we can envisage a number of improvements. 
We have assumed the purely radial expansion (with a given velocity R{t) of the plasma boundary) and hence the shape of the 
plasma (in the xy-plane) remains isotropic during plasma expansion. For realistic laser-generated plasmas this is valid when 
the kinetic energy density of the expanding plasma much exceeds the magnetic field energy density |[ll|2|l- However, after some 
period of accelerated motion, the plasma gets decelerated as a result of the external Lorentz force acting inward and the above 
mentioned condition may be violated at the later stages (t ^ td) of the expansion. In this case the Lorentz force density which 
is anisotropic in general, cannot be neglected in Eq. ([T]i. Thus one can expect some deformation of the initially isotropic plasma 
surface and the plasma radius R{t, (p) should be treated as a function of (p. On the other hand, in this regime the plasma 
radius R{t, (p) cannot be treated as a given function of time and should be determined self-consistently using, for instance, the 
equation of balance of plasma (^ pR^) and magnetic field (~ H^) energies |T9t]. Thus, plasma dynamics can be described (at 
least qualitatively) inserting Eqs. (l33Tl-(l35]l and (l38T l. (|39] l into the energy balance equation which yields a first-order differential 
equation for R{t, ip). Within the scope of this analysis, the deformation of the plasma surface as well as the initial stage of 
plasma acceleration, the later stage of deceleration and the process of stopping at the point of maximum expansion could be 
examined numerically. We intend to address these issues in our forthcoming investigations. 
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Appendix A: Some summation formulas involving the Bessel functions 



Using the Fourier-Bessel expansion of an arbitrary function of a real variable 112511 one can derive some summation formulas 
involving the Bessel functions which are used in the main text of the paper. The first relation is obtained by considering the 
Fourier-Bessel series for the quadratic function in the interval ^ x < 1, 

where A„ (with rt = 1, 2, ...) are the positive zeros of the Bessel function, Jo(A„) = 0, arranged in ascending order of magnitude 

il. 

On the other hand using the known summation formula [25] (valid at ^ x < 1) 



one can represent Eq. (lAlt in another form 



V -^"^^"^^ = i rA2^ 
^A„Ji(A„) 2 ^ ' 



This latter relation is valid at ^ x ^ 1. 

Next taking the x-derivatives of Eqs. dAlb and (IA2b one obtains 

respectively. The latter formula is valid at ^ a; < 1. Therefore, using the second relation in Eq. (IA4) the first summation 
formula in (IA4b can be represented in the form 
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Finally, substituting a; = 1 in Eq. iA5\ we arrive at 



A2 ~ 4 ' 

n=l " 



(A6) 



Next we derive another kind of summation formulas involving the zeros of the Bessel function Jo{z). For that purpose 
consider the known summation formula 1,25.1 



lojxz) _ „ A„Jo(Ana::) 



(A7) 



where In{z) (with n ~ 1,2, ...) is the modified Bessel function of the first kind, z and x (with ^ a; < 1) are real variables. 
Differentiating Eq. (IA7b with respect to x and using the second relation in Eq. ( IA4l i we arrive at 



E 



Ji(A„x) 



hjxz) 

2zIo{z) ^^{Xl+z2)J,^Xn) 



which, in particular, at x = 1 yields 



CXD ^ 

E T2~T 



h{z) 

-A2+z2 2zIo{z)' 



(A8) 



(A9) 



This relation can be used for deriving other summation formulas. For instance, extracting Eq. ( |A9l l from Eq. (IA6) and after some 
manipulations one obtains 



-^A2(A2z + l) 
4 



,^^A2(A2z + l) 



2 - 



2/i(C) 

ao(c) ' 

4/i(C) 

ao(c) 



^i(C) 



4(C) 



(AlO) 



(All) 



Here ^ = 1/ y/z- We also mention the asymptotic behavior of the functions Si (z) and .^2 (z). At small argument, z <C 1, these 
functions behave as Si(z) ~ '^2{z) — 1 while at z 1 they decay as Si(z) ~ l/8z, S2(z) ~ 1/48^^. 
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